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ABSTRACT 

The recent detections of TeV gamma-rays from compact binary systems show that 
relativistic outflows (jets or winds) are sites of effective acceleration of particles up to 
multi-TeV energies. In this paper, we discuss the conditions of acceleration and radia- 
tion of ultra-relativistic electrons in LS 5039, the gamma-ray emitting binary system 
for which the highest quality TeV data are available. Assuming that the gamma-ray 
emitter is a jet-like structure, we performed detailed numerical calculations of the en- 
ergy spectrum and lightcurves accounting for the acceleration efficiency, the location 
of the accelerator, the speed of the emitting flow, the inclination angle of the system, 
as well as specific features related to anisotropic inverse Compton (IC) scattering and 
pair production. We conclude that the accelerator should not be deep inside the binary 
system unless we assume a very efficient acceleration rate. We show that within the IC 
scenario both the gamma-ray spectrum and flux are strongly orbital phase dependent. 
Formally, our model can reproduce, for specific sets of parameter values, the energy 
spectrum of gamma-rays reported by HESS for wide orbital phase intervals. However, 
the physical properties of the source can be constrained only by observations capable 
of providing detailed energy spectra for narrow orbital phase intervals (A0 <^ 0.1). 



1 INTRODUCTION 



sion of LS 5039 based on both leptonic and hadronic inter- 



LS 5039 is a binary system consisting of a very bright star 
and a compact object (a neutron star/pulsar or a black- 
hole). The source emits X-rays (see e.g. |Bosch-Ramon| 



actions (e.g. Dermer & Bottcher (20061; Aharonian et al 



(2006b I; Paredes et al. (20061; Bednarek (2007l). At the 



et al. (20071 and references therein; Goldoni et al. (20061; 



De Rosa et al. (20061), and presumably also MeV (Strong 
et al.||2001| )~and GeV ( [Paredes et"al][2000l [2002] ) gamma- 
rays. Recently, LS 5039 has been detected in very high en- 
ergy (VHE) gamma-rays by the HESS array of atmospheric 
Cherenkov telescopes ( [Aharonian et alT]|2005[ ). The TeV ra- 
diation of the source is clearly modulated with a period 
3.9078 ± 0.0015 days ( [Aharonian et al. [ [2006a!), w hich per- 
fectly coincides with the orbital period of the object ( [Gas ares] 
et al.|2005| ). 



same time, a scenario in which the nonthermal radiation is 
related to an ultrarelativistic pulsar wind, remains an al- 
ternative option (Dubus 2006a). In this scenario, LS 5039 
would behave in a similar manner to the binary pulsar sys- 
tems (e.g. PSR B1259-63/SS2883), where the production 
site of gamma-rays is (most likely) related to the pulsar 



wind termination shock (Maraschi & Treves 1981 



& Arons|[1997[ [Kirk et al.[[1999' 'Khangulyan et al 



Tavani 



20071 



The nature of the compact object in LS 5039 is not 
yet firmly established, given the uncertainty related to the 
orbital inclination angle, 13° < i < 64°, which does not al- 
low precise estimate of its mass. The optical line analysis, 
together with the assumption of orbital pseudosynchroniza- 
tion favour a rather small inclination angle, i ~ 25° ( [Gas ares] 
et al. 20051. This would be an indication of a rather high 
mass for the compact object, M ~ 3.7 Mq, and therefore 
could be interpreted as an evidence for its black-hole na- 
ture. In such a case, one would expect the realisation of 
the microquasar scenario, in which the nonthermal processes 
take place in a jet related to an accreting compact object. 
The original classification of LS 5039 as a microquasar was 



Whereas in the microquasar scenario, the particle acceler- 
ation and gamma-ray production are possible throughout 
the entire jet (i.e. both inside and outside of the binary 
system), in the standard pulsar wind model the particle 
acceleration and radiation take place well inside the system 
( Dubus [[2006a I, namely at distances from the compact ob- 
ject much smaller than the s eparation between th e stars, i.e. 
Rorh = (1-4 - 2.9) • lO^^cm ( [Gasares et al.|2005 l. 



based on its extended radio emission features ( Paredes et al, 
|2000|). In the framework of this scenario, a number of mod- 



els have been proposed to explain the TeV gamma-ray emis- 



The discovery of modulated TeV gamma-ray emission 
in LS 5039 is a strong indication that gamma-rays are pro- 
duced close to the binary system. On the other hand, if the 
gamma-ray production region is located inside or very close 
to the system, we expect distinct photon-photon absorption 
features caused by the interaction between the gamma-rays 
and the stellar radiation field. Within the leptonic models 
for gamma-ray production, two other effects have a strong, 
direct or indirect, impact on the formation of the gamma-ray 
spectrum. One effect is related to the anisotropic IG scat- 
tering; the other effect is related to the maximum energy of 
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Figure 1. Sketch of the system. 

relativistic electrons determined by the balance between the 
acceleration and radiative cooling rates. Finally, the produc- 
tion and gamma-ray absorption processes should be coupled 
with an appropriate treatment of the particle cooling and 
propagation (e.g. particle diffusion and advection along the 
jet). Thus, in the context of leptonic models, the spectral 
shape of the gamma-ray emission, in particular its orbital 
phase-dependence, contains unique information about the 
location of the particle acceleration and gamma-ray produc- 
tion regions. 

In what follows, we show that the spectral and temporal 
features of the VHE gamma-radiation reported by HESS 
from LS 5039 can be reproduced with a leptonic model that 
invokes anisotropic IC scattering of electrons and absorption 
of gamma-rays under certain requirements concerning the 
location and efficiency of the electron accelerator. 



2 PHYSICAL PROCESSES IN THE SYSTEM 

2.1 Radiation and acceleration processes 

The gamma-ray spectrum of LS 5039 extends to very high 
energies, up to ~ 10 TeV and beyond ( [Aharonian et alT| 
2006a I. Since the parent electrons should have even higher 
the acceleration of > 10 TeV electrons requires 



energies, 

rather special conditions, especially if the radiation is formed 
inside the binary system, where electrons suffer strong IC 



energy losses (Aharonian et al. 2006b I. Assuming for sim- 
plicity that the gamma-ray emitter is a jet-like structurcj^ 
perpendicular to the orbital plane, one derives the energy 
density of the stellar radiation field in the emitting region: 

^ The assumption on the jet-Uke structure here does not nec- 
essarily mean adopting the microquasar scenario (i.e. an accret- 
ing system with a jet). In fact, a jet-like structure could also be 
produced in other types of object, e.g. in binary pulsar systems 
i Bogovalov et al.|2007 l. 




B, G 

Figure 2. Contour plot of the maximum energy of accelerated 
electrons in the B-d plane for different values of the rj parameter: 
Tj = 1 (solid lines), rj = 10 (dashed lines), = 100 (dotted lines). 
Thick and thin lines correspond to 30 TeV and 1 TeV electrons, 
respectively. 



Wr — Li, / A-k{R^^^ + Z^)c, where L* is the star luminosity 
and Z is the distance to the compact object. 

A sketch of the geometry of the system is shown in 
Fig. [l] The radiation of the companion star in LS 5039 
is characterised by a black-body spectrum of tem perature 
kT « 3.3 eV and luminosity « 7 • 10'^* erg/s (Casares 



et al. 2005 1 . Since Z cannot be much larger than the or- 



bital distance (otherwise the periodic component of radia- 
tion would be smeared out), the energy density of the target 
photons typically varies between 10 erg/cm^ at Z ~ 10^'^ cm, 
and 1000 erg/cm'^ at the base of the jet around periastron. 

The IC scattering of TeV electrons on the starlight pho- 
tons of average energy 3kT « 10 eV takes place deep in the 
Klein-Nishina (KN) regime. The characteristic cooling time 
of electrons in this regime can be approximated with good 
accuracy as ( Aharonian et al.|2006b l: 

tKN ~ L7 • lO'wcT'S" TeV S 



(1) 



where Eo tcV = E/lTeV is the electron energy in TeV units, 
Wo = Wi /lOO erg/cm^ and di3 = V-R^t + ^VlO" 

cm. The 

maximum acceleration energy of electrons is achieved when 
the cooling time approaches the acceleration time (face), 
which it is convenient to present in the following general 
form: 



: ryr-L/c ^ 0.1 Ec tcV Bq q 77 s , 



(2) 



where tl = E/eBo is the Larmor radius, and Bo G = 
Bo /I G is the strength of the magnetic field in the ac- 
celerator. The parameter rj characterises the acceleration 
efficiency. Generally, rj ^ 1, and only in so-called extreme 
accelerators r/ ^ 1 ( Aharonian et aL]|2002 1 . From the con- 
dition tacc = tKN one obtains: 



4 • 10'" [Bo GV' 



TeV. 



(3) 



Formally, for an extreme accelerator [i] < 10) and Bq ^ 
0.3 G, the maximum energy of electrons can, under KN IC 
energy losses, exceed 100 TeV even for a radiation energy 
density Wr — 10^ erg/cm^. However, for such a magnetic 
field strength in the accelerator, synchrotron energy losses 
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Figure 3. The same as in Fig.[2]but only for r} = 10. 

dominate over the Compton losses. The synchrotron cooling 
time is: 



lO^B, 



2 p-l 
G-^G TcV 



thus the condition face ~ tsy gives: 



-1/2, 
-"0 G 



(4) 



(5) 



Another fundamental condition, tl < Zo, where Zq is 
the location of the accelerator and taken as an upper-limit 
of its Hnear size (ia 12 ~ gives: 



Eo < 3-10" /a 12-Bo G TeV 



Ee <3-lo'2'o/i2BoGTeV,(6) 



where Zo/12 = Zo/10^^. 

In Fig. [2] we show the maximum-energy contour plot in 
the parameter plane {d, B) based on the conditions given 
by Eqs. (|3j[5][6|. The curves correspond to three values of 
the r] parameter: 77 = 1, 10 and 100. It is seen that elec- 
trons can be accelerated up to ~ 30 TeV energies deep in- 
side the binary system (d ~ 2 ■ lO^'^ cm, i.e. Zq <^ Roih) 
only in the case of an extreme accelerator with r) < 10. It is 
worth noting that in dense radiation environments the accel- 
eration efficiency can be significantly enhanced through the 



so-called converter mechanism (Derishev et al. 2003 Stern 
2003). Therefore, quite small values of the acceleration effi- 
ciency, i.e. 77 ~ 1, cannot be a priori excluded. 

In Fig. |3j we show the maximum electron energy map 
in the (d, B) plane for rj = 10. It is seen that for r\ > 10, 
electrons can be accelerated to energies > 30 TeV only in 
an environment with Bo ^ 0.3 G located at Zo > 10^^ cm. 

The hardness of the HESS reported spectra of gamma- 
rays allows us to put constrains on the magnetic field 
strength in the emitter. Namely, the reported photon in- 
dices, ranging from ~ 2 — 2.5, indicate that the magnetic 
field energy density in the emitter should be significantly 
smaller than the target photon energy density. Indeed, if 
these gamma-rays are produced via IC scattering, the elec- 
tron differential spectrum must be harder than oc E^'^. For 
dominant synchrotron cooling, the electron energy distribu- 
tion at high energies is softer than oc E~^. Even in case 
of a monoenergetic injected electron spectrum, synchrotron 
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Figure 4. Electron propagation length versus location of injec- 
tion in the system under different magnetic fields. Calculations are 
performed for a black-body radiation field with T = 3.8 ■ 10* K, 
with a radiation energy density corresponding to that of the star 
at distance d (i.e. dilution coefficient g = (H*/2d)^ with = 
6.5 ■ 10^^ cm), and for three values of the magnetic field B = 0.03, 
0.1, 0.3 G, which are constant along the propagation path. The 
advection velocity was assumed to be Vadv = 10^" cm/s. We 
note that low energy electrons (100 MeV) can propagate very far, 
whereas 100 GeV electrons have the shortest propagation length 
and, depending on B, 10 TeV electrons can propagate farther 
than 1 TeV electrons. 

cooling results in a oc E^^ type electron energy distribution. 
However, IC energy losses taking place in the KN regime al- 
low such a harder than oc _E~^ electron energy distribution. 
Thus, hard VHE spectra require Ikn < tsy, or 

Bg < o.ed^iE;'^^^ (7) 

for the magnetic field in the emitter. 
2.2 Propagation of electrons 

Relativistic electrons can propagate along the jet. In such 
a case, their energy cooling proceeds under changing phys- 
ical conditions. This can have a strong impact on the re- 
sulting gamma-ray radiation. We consider here two possible 
transport mechanisms: diffusion and advection along the jet. 
The propagation length depends on the diffusion coefficient 
D, on the bulk velocity of the jet (or advection velocity) 
Vadv, and on the radiation cooling time tcooi. The particle 
transport is important when the propagation distance AZ 
is comparable to the separation between the injection point 
and the compact object (Zo). Thus, the particle transport 
can be described by a dimensionless parameter k = AZ/Zq, 
which is dominated, depending on the propagation regime, 
by KdiEF = Adiff/.^o or Kadv — Vadvicooi/-^o. The diffusion 
length Adiff = ^/273^^, thus: 



Kdiff- = 0.1 Zq^-^^^E, 



1/2 



TcV^G 



1/2 ft 
VIO2 



^s^ V 



D 



1/2 



(8) 



where Deohm is the diffusion coefficient in the Bohm regime, 
and Eb TeV is the electron energy in TeV units. Thus, elec- 
tron diffusion has no impact on the TeV radiation, unless 
the diffusion is far from the Bohm regime {D 3> -DBohm) or 
the magnetic field is very small in the emitter. For advection 
we obtain: 
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Figure 5. Dependence of the IC energy spectrum on the interac- 
tion angle: die = O.Itt, 0.27r, O.Stt, O.Stt, O.Ttt, tt. It is assumed that 
monoenergetic electrons (_Ec = 0.1 TeV) interact with black-body 
photons (T = 3.8 ■ lO" K). 



10 



10"'cm/s 



1 J 

0/12-' 



(9) 



when energy losses are dominated by IC scattering, and 



10iOcm/s^o/i2^"TcV^G ^ 



(10) 



for dominant synchrotron energy losses. Therefore, for a 
mildly relativistic outflow (\4dv ~ lO'^" cm/s), k can eas- 
ily exceed 1 under dominant KN IC energy losses. In Fig.|4j 
the results of the electron advection length calculations are 
shown for four different energies = 100 MeV, 100 GeV, 
1 TeV, and 10 TeV. Here we consider a Vkdv = 10^° cm/s. 
Three constant values of the magnetic field were considered: 
B = 0.03, 0.1, and 0.3 G. The radiation field was assumed 
to be black-body with T — 3.8 • 10* K, and energy density 
corresponding to that of the star at distance d (i.e. dilution 
coefficient g — (i?*/2d)^, where _R* = 6.5- 10^^ cm). As seen 
in Fig. |4j electrons with energy Ec ~ 100 GeV have the 
shortest advection length. Due to cooling in the Thomson 
regime (tx oc E^^), electrons with smaller energy have a 
longer propagation length. In case of higher electron energy 
{Ec > 100 GeV) , IC cooling proceeds in KN regime, leading 
to an increase of the propagation length with the growth of 
the electron energy unless synchrotron energy losses start to 
dominate. All this is clearly seen in Fig. [4] 

As discussed above, the parameter k varies over a 
rather wide range depending on the electron energy and 
the accelerator-companion star separation distance. Thus, 
high energy electrons can propagate far away from the ac- 
celeration site. In the next section we show that this effect 
can significantly change the gamma-gamma optical depth 
for the produced VHE radiation and the IC scattering angle 
between electrons and target photons. 



2.3 IC and gamma-gamma absorption 

Because of the strong anisotropy and inhomogeneity of 
the target photon distribution, IC scattering and gamma- 
gamma absorption depend on the electron location and mo- 
mentum direction. Regarding IC scattering, the interaction 
angle die has a strong impact on the formation of radia- 
tion as shown in the past by several authors in different as- 



trophysical scenarios (e.g. 


Bogovalov & Aharonian (20001; 


Khangulyan & Aharonian 


(|2005a|b|); |Dermer & Bottcher 



(2006 



bution (SED) produced by a monoenergetic distribution of 
electrons of i?e = 0.1 TeV interacting with photons with a 
black-body distribution (T = 3.8 • 10* K) at different an- 
gles 6^10 ranging from O.lvr to tt. It is seen that, depending 
on 6*10, for electron energies Ec > m^c^/kT, the IC scat- 
tering can proceed either in the KN {E-,/Ec ~ 1) or in the 
Thomson [E-y/Ee <^ 1) regimes. In In Fig. |6] we show the 
computed SED corresponding to the IC scattering between 
electrons with a powerlaw energy distribution and photons 
with a black-body distribution (T = 3.8 • 10* K). Results 
are shown for different 9ic ranging from O.Itt to tt. Since the 
interaction angle strongly varies along the orbit and with 
Z, the final gamma-ray spectra significantly depend on the 
orbital phase and location of the emitter. Fig. [7] we show a 
two-dimensional representation of the IC scattering proba- 
bility in the plane (7,6/7; where 7 is electron Lorentz fac- 
tor and e is the outgoing photon energy) for three different 
9ic ~ 5°, 90° and 175°, which are shown altogether with the 
average cross-section case. Depending on the electron distri- 
bution function, the IC anisotropy could have a different 
impact on the gamma-ray spectrum. In the case of a mo- 
noenergetic or narrow electron distribution (e.g. (5-function) 
around 0.1 TeV (~ m^c^/kT) the IC scattering at high in- 
teraction angles results in higher fluxes and harder spectral 
shapes (see Fig. [5|. For broad (e.g. power-law) electron dis- 
tributions, IC scattering at small interaction angles leads to 
a rather low flux, but the spectral shape is harder than for 
large scattering angles. The hard spectral shape is due to a 
transition from the Thomson, in which the scattering is sen- 
sitive to the interaction angle, to the KN regime, in which 
the angular dependence is weak. We note that the smaller 
the interaction angle, the larger the transition energy. 

The gamma-ray absorption in close binaries was studied 
by, e.g ., [Moskalenko fc Karaku"Ia| ( |1994[ ), |B6ttcher fc Dermer| 



(20051 and Dubus (2006b I. In general, the gamma-gamma 
optical depth depends on the distance d between the emit- 
ting electrons and the companion star, and on the scattering 
angle (assuming a point-like source of target photons) Oic'- 



(11) 



where 6 (I) is the angle between the directions of the 

gamma-ray and the stellar photon at distance / from the 
emitting point. The absorption probability strongly depends 
on the interaction angle. This is shown in Fig. |8] where the 
probability of gamma-gamma absorption on a monoener- 
getic distribution of target photons (e = 10 eV) is plotted 
for the interaction angle = O.Itt, 0.2tt, 0.3tt, 0.5tt, 0.7tt, tt. 



T{d,6ic,E^)= J dl {1 — cos 9^.y) 

line of sight 

X / de app{eE^{l ~ cos6^-,))nph{e). 
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Figure 6. Angular dependence of the IC cross-section. The re- 
sult of the interaction of a powerlaw electron distribution with 
a black-body distribution of target photons (T = 3.8 ■ 10'' K) 
is shown for the following values of the interaction angle: diQ = 
O.Itt, 0.27r, O.Stt, O.Stt. 0.77r, tt. 
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Figure 8. The probability of gamma-gamma absorption for 
a monoenergetic distribution of target photons (e = 10 eV) 
for the following values of the interaction angle: 9-^^ = 
O.Itt, 0.2-k, O.Stt, O.Stt. O.Ttt, tt. 
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Figure 7. A 2D representation in the (7,e/7) plane of the IC 
interaction probability for different = 5° (upper-left panel), 90° 
(upper-right panel), 175° (lower-left panel). The same is shown for 
the angle averaged IC interaction probability (lower-right panel). 



Concerning the d-dependence of the optical depth, in the 
point-like approximation for the geometry of the target pho- 
ton source, the gamma-gamma optical depths for arbitrary 
do and d can be scaled as: 



-jT{do, 9ic, E-, 



(12) 



From previous considerations, the smaller 6^10 the 
harder the produced radiation. In addition, since 9ic ^ 1 
implies also 6.^.^ <C 1 in Eq.(lll, the threshold energy 



for gamma-gamma absorption significantly increases up to 
^-th 3> 100 GeV. Both effects result in hard and unabsorbed 
gamma-ray spectrum above 1 TeV. 



2.4 Radiation of secondary electrons 

The interaction of gamma-rays with the stellar radiation 
field produces a population of secondary electron/positron 
pairs (secondary electrons hereafter) inside the binary sys- 
tem. A significant fraction of the primary gamma-ray energy 
can be transfered to these particles. In Fig. [9] we show an 
example of the fraction of gamma-ray energy absorbed in 
the radiation field of the star. Assuming an isotropic point- 
like emitter located at distance d from the optical star, we 
calculate the absorbed fraction of primary gamma-ray en- 
ergy in the energy range 0.1 — 10 TeV (for a E^'^ spectrum 
of primaries). One can see that, even for an emitter located 
rather far from the star, radiation of secondaries can signifi- 
cantly contribute to the observed emission from the source. 
The secondary electrons radiate via the synchrotron or IC 
channels. The ratio of these fluxes depends on the energy 
band, the magnetic fleld strength and the distance to the 
companion star. This effect may lead to the formation of an 
electromagnetic cascade as long as the energy of electrons 
and gamma-rays exceeds a certain energy: 



E > 



GeV 



(13) 



For the effective development of the cascade, IC energy 
losses should dominate over synchrotron energy losses. From 
Eqs. ([T| and (|4| one obtains: 

T,V < 2 • IQ-^WrB-^ , (14) 



E, < 1 



TeV 



(15) 



where wb = B'^/8tv is the energy density of the ambient 
magnetic fleld in the system, which is distinguished here 
from the magnetic fleld inside the emitter. This relation al- 
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Figure 9. Fraction of the primary gamma-ray energy emitted in 
the range 0.1 - 10 TeV (-Eaba) which was absorbed as a function 
of the distance to the star in stellar radius units. We assumed a 
power-low distribution of primary gamma-rays (oc E-^). 



lows us to introduce a critical value of the magnetic field in 
the binary system: 



Ba c = 10 



is 



7. 1038£££ 



1/2 



(16) 



where R is the distance between the gamma-ray absorption 
point and the star. If the magnetic field is stronger than 
_Ba c, the electromagnetic cascade will not affect the energy 
range i5o > 1 TeV, implying a purely absorbed spectrum for 
> 1 TeV. 

Although the magnetic field within the system is un- 
known, the magnetic field at the surface of O type stars 
can be typically ~ 100 G - 1 kG (see e.g. Usov & Melrose 



(19921; Donati et al. (2002l). Taking into account the value 
of Ba, c (see Eq. (16l), efficient cascading at > 1 TeV seems 
unlikely. For such magnetic fields, synchrotron emission is 
the main radiation channel at these energies, and electro- 
magnetic cascading will be suppressed. 



3 MODEL DESCRIPTION 

In this work, we adopt a model for gamma-ray emission in 
LS 5039 in which emission is produced in a jet-like structure 
perpendicular to the orbital plane. Also, we assume that the 
acceleration takes place in a compact region located in the 
jet at Zq ( although the emitting region could be much larger 
because of particle advection) and Zq does not depend on 
the orbital phase. In Fig. [l] the accelerator and the emitter 
are shown as a black and a grey region, respectively. 

The acceleration/cooling time-scales in the system are 
much shorter than the orbital period. Therefore, the electron 
distribution function n can be described by a steady-state 
equation. In the case of an one-dimensional jet, n(7, Z) is 
determined by the following equation: 



Vadv 



dn 97(7, Z)n 



10' 



10' 
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Figure 10. Calculated SED for three different Bq = 0.03, 0.1 
and 1 G. The rest of parameters have been fixed to the next 
values: Zq = 2 ■ lO^^ cm, t? = 1, V^dv = 10^° cm/s, and i = 25°. 
Gamma-ray absorption is included. 



where the Vkdv is the jet bulk velocity; 7(7, Z) is the elec- 
tron energy loss rate, which accounts for IC and synchrotron 
mechanisms and depends both on the electron energy and 
position in the jet; Zq is the location of the accelerator; and 
5(7) is the injection electron spectrum. We note that 5(7) 
and 7 are slow functions of time, i.e. they vary adiabatically 
along the orbit. Because of a lack of knowledge concerning 
the acceleration processes occurring in LS 5039, we adopt 
a phenomenological function for 5(7). Namely, we assume a 
power-law function with an exponential high-energy cutoff: 



9(7) = Qo 7 " exp{-7mc^/£o} 



(18) 



-Eo is the electron cutoff energy determined from the balance 
between acceleration and energy loss rates (see Sec. |2.1[ ); a 
is fixed to 2, which is a common value for injected spectra 
of non-thermal electrons, and Qo is a constant. 

We have numerically solved Eq. (17 1 taking into ac- 



dz 



97 



V.d.q{'y)S{Z - Zo) 



(17) 



count IC and synchrotron radiation processes. We neglect 
at this stage ionization or Bremsstrahlung losses because 
they are not relevant in the present context. We assume 
that the main targets for IC scattering are provided by the 
star, and that the magnetic field decreases with the dis- 
tance as B = Bo{Z/Zo)~^ . In order to accelerate electrons 
to energies >10 TeV the magnetic field should be within 
-Bo ~ 0.01 — 1 G, depending on the acceleration rate and 
location (see Fig. [2|. 

To study the impact of the magnetic field on the IC 
and synchrotron radiation components, three different Bo ~ 
0.03, 0.1 and 1 G have been adopted. The rest of relevant 
parameters have been fixed to the values: Zq = 2 • 10^^ cm, 
1] = 1, Vadv = 10^" cm/s, and i — 25°. The correspond- 
ing broadband SED is shown in Fig. |10| We note a very 
hard synchrotron spectra in the interval 1 keV - 1 MeV 
(with a spectral energy index ~ 0.5). This is explained 
by the hardening of electron energy distribution due to 
the radiative cooling in KN regime ( [Khangulyan fc Aharo-| 
nian||20()5a Moderski et al.|2005l. Since for magnetic fields 
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Figure 11. The distribution of electrons along the jet for different 
energies. B = 0.2{Z/Zo)-^ G, Zq = 2 ■ lO^^ cm, and Vadv = 
fO^O cm/s. 



B < 1 G the energy losses are dominated by IC scattering, 
the synchrotron radiation fiux is significantly less than the 
IC gamma-ray flux, and below also the observed fluxes at 
X-rays ( [Bosch-Ramon et al.]|2007 l. 

It is worth discussing the importance of the electron 
propagation, already commented on in Sec. |2.2| Namely, we 
expect the most efficient propagation of electrons at very low 
(Thomson regime) and at very high (deep KN regime) ener- 
gies as long as IC energy losses dominate over synchrotron 
ones. In Fig. we show the distribution of electrons along 
the jet. It is seen that, the higher the energy, the larger the 
propagation length of electrons. This is simply the result of 
the KN effect. However, for the assumed magnetic field, this 
tendency terminates around 30 TeV, when the IC energy 
losses are overcome by the synchrotron ones. We note also 
that, below Ec < 10 GeV, the lower the energy, the longer 
the propagation length of electrons, because in this energy 
interval the IC scattering occurs in the Thompson regime. 

In fact, electrons at large Z will interact with weaker 
magnetic and stellar radiation fields, thus it is convenient 
to introduce a parameter which characterises the effective 
electron average location: 



J {Zn{Z,E^)/if)dZ 

Za 

OO 

/ {n{Z,E^)/d^)AZ 



(19) 



The .Zo-dependence of thi s effective electron average loca- 
tion is shown in Fig. 12 where the l/d^-weighted aver- 



age locations of the electrons of different energy are shown 



for three different = 1 ■ 10^^, 2 
B = 0.4(Z/10^^cml-i G and 14dv 



10^% and 5 • 10'^ cm. 
10^" cm/s are adopted 
for all three casea^ The 1/d^-weight comes from the fact 
that the stellar radiation field decreases as 1/d^, giving an 
idea on where electrons of a certain energy will radiate the 



We recall (see Eq. [8| that diffusion is negligible for these pa- 
rameter values. 




100 



E, TeV 

Figure 12. The electron -weighted average location for three 
different Zq, i.e. 1- lO^^^ 2-10^'^, and S-IO^^ cm, adopting V^dv = 
IQiO cm/s and B = OA{z/lO^'^cm)-^ G. 



most via IC. We see that Z{Ec) significantly depends not 
only on the energy but also on the location of the accelera- 
tor. 

This rather unusual energy-dependence of the propaga- 
tion of electrons in the jet can be reflected in the features 
of the electron radiation. We note, however, that the spa- 
tial distribution of the gamma-rays along the jet does not 
simply repeat the distribution of electrons. Two additional 
effects, related to the interaction angle in the IC scattering 
and gamma-gamma absorption, lead to a significant modi- 
fication of the gamma-ray effective distribution along the 
jet. These effects are shown in Figs. |13| and |14[ respectively. 

Finally, we want to emphasise the impact of the angular 
dependence of the IC scattering on the observed spectrum. 
The importance of this effect for the IC SED is demonstrated 
by Fig. |15| The two curves shown in this figure are obtained 
under different assumptions. Whereas curve 1 is calculated 
using the angle-averaged Compton cross-section, curve 2 is 
calculated using the precise angular dependent cross-section 
and the exact interaction geometry, being averaged over the 
orbital phase interval 0.45 < < 0.9 (see Sec.[4]|. 



4 RESULTS 

We now apply the model developed in the previous section 
to LS 5039. We have computed the SED and the orbital 
lightcurves for different relevant situations. Some of the 
SED presented here are calculated for the inferior conjunc- 
tion of the compact object, i.e. — 0.72 (infc; when the 
compact object is between the observer and the compan- 
ion) and the superior conjunction, i.e. — 0.06 (supc; when 
the compact object is behind the companion star with re- 
spect to the observer). HESS has reported on the temporal 
and spectral characteristics of the emission from LS 5039 
(jAharonian et al. 2006a). Because of the lack of statistics, 
HESS has reported the SED for two different phase inter- 
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Figure 13. The averaged location for the production of 10 TeV 
gamma-rays produced in the jet. Zq = 2 - 10^^ cm and i = 25° are 
adopted. This shows the importance of the IC interaction angle 
for the resulting gamma-ray production. 



vals: 0.45 < < 0.9 (~infc); tj) < 0.45 and (p > 0.9 (~supc). 
In order to explore the model parameter space by comparing 
the calculations with the HESS data, some of the calculated 
SED are averaged over the same phase intervals. 



4.1 Exploring the model parameter space 

In this section, we fix Bo — 0.05 G, which is weak enough 
to provide a radiation dominated environment for ~ TeV 
electrons (see Eq.([7|). We also assume a high acceleration 
efficiency with rj = 10, required by the extension of the HESS 
reported spectrum well beyond 10 TeV (see Fig. [2|. Given 
the discussions in previous sections, we assume the following 
values for the remaining parameters: Zo = 10^^, 2 ■ 10^^, 5 • 
10^2 cm; Kdv = 10^ 10^" cm/s; and i = 15°, 25°, 55°. We 
compute the gamma-ray SED for infc and supc. These are 
two distinct phases with very different physical conditions 
for gamma-ray production and absorption. For instance, the 
maximum gamma-ray absorption (and the minimum Eth of 
this process) occurs in supc. 

In Fig. |16| we show the gamma-ray SED when the elec- 
tron propagation is negligible, i.e. 14dv = 10^ cm/s, for two 
different i = 15°, 55°, and Zq = 2 ■ 10^^ cm. This relatively 
slow advection velocity implies that the gamma-rays are pro- 
duced at Zo- As seen from Fig. |16| for small inclination an- 
gles (left panels), there are no strong spectral differences 
between supc and infc, and the orbital lightcurve is rather 
flat, as shown in Fig. |17| (left panel). This is caused by the 
small changes in the angle of interaction for i = 15°. We 
note that the orbital lightcurve would peak more strongly 
around apastron passage for slightly larger -q, because of the 
strong r;-dependence of the highest energy gamma-rays for 
the adopted Zq. We note that the counter-jet impact is not 
strong for low inclination angles but below ~ 10 GeV, be- 
cause its VHE emission is strongly absorbed. This jet dom- 
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Figure 14. The gamma-gamma absorption factor, exp (— r), cor- 
responding to three different energies for the gamma-ray produced 
in the jet: 0.1 (solid lines) 1 (dashed lines), and 10 TeV (dotted 
lines). Thin lines represent the absorption factor of gamma-rays 
produced in the accelerator, at Zq = 2 ■ 10^^ cm. The absorption 
factors relevant for the average location of the gamma-ray emit- 
ter are shown with thick lines. The inclination angle is assumed 
to be i = 25°. This shows the importance of the gamma-gamma 
interaction angle for the gamma-ray absorption. 
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Figure 15. The IC SED obtained with two different IC cross- 
sections. Curve 1 is calculated using the angle averaged inverse 
Compton cross-section. Curve 2 corresponds to calculations in 
which the angle-dependent cross-section has been averaged only 
over the orbital phase interval 0.45 < </> < 0.9, roughly when the 
emitter is between the companion star and the observer (i.e. 8ic 
relatively small; see Sec.[4]l. The calculations were performed for 
Zo = 2 ■ 10^2 B = 0.1(Z/Zo)-l G, r; = 35, i = 25°, and 
Vadv = 10^" cm/s. We show added both the jet and the counter- 
jet component. 
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Figure 16. Gamma-ray SED for supc (bottom panels) and infc 
(top panels) for different inclination angles: i = 15° (left panels) 
and 55° (right panels). The remaining parameter values are Zq = 
2 • 10l2 cm, y^dv = 10^ cm/s, = 10 and B = 0.05(Z/Zo)-^ G. 
The jet and the counter-jet components, plus the summation of 
both, are shown. The absorbed components are presented with 
solid lines, and the production components are shown with dotted 
lines. 




Figure 17. Gamma-ray orbital lightcurves for the same cases 
and parameter values as in Fig. |16| The jet and the counter-jet 
components, plus the summation of both, are shown. The ab- 
sorbed components are presented with solid lines, and the pro- 
duction components are shown with dotted lines. We note that 
the counter-jet production component is not seen since it is off 
scale. 



inance is a purely geometrical effect in the IC interaction 
and gamma-gamma absorption. For high inclination angles, 
the interaction and absorption angles change significantly 
yielding strong spectral and flux variations along the orbit, 



as shown in Fig. 17 (right panel). Two remarkable features 
are a double-bump orbital lightcurve, which is already sig- 
nificant at i « 25°, and a more relevant VHE counter-jet 
for i > 25° due to weaker gamma-ray absorption for this 
component. 

For an advection velocity 14dv = 10^ cm/s and a fixed 
i = 25°, the dependence of the results on Zq is demonstrated 
"■^^ cm (left panels), fast IC energy 



by Fig. 



18 



For Zo = W 

losses around periastron prevent efficient electron accelera- 
tion beyond 1 TeV unless < 10. A result of this effect is 
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Figure 18. Gamma-ray SED for supc (bottom panels) and infc 
(top panels) for different Zq = 10^^ cm (left panels) and 5 ■ 10^'^ 
(right panels). The rest of parameters are assumed to be i = 25°, 
Vadv = 10^ cm/s, rj = 10, and B = 0.05(Z/Zo)"^ G. The jet 
and the counter-jet components, plus the summation of both, are 
shown. The absorbed components are presented with solid lines, 
and the production components are shown with dotted lines. 




Figure 19. Gamma-ray orbital lightcurves for the same cases 
and parameter values as in Fig. |18| The jet and the counter-jet 
components, plus the summation of both, are shown. The ab- 
sorbed components are presented with solid lines, and the pro- 
duction components are shown with dotted lines. We note that 



the counter-jet production component for Zq ■ 
seen since it is off scale. 



: 5 • 10 cm is not 



a strong steepening in the highest energy part of the unab- 
sorbed SED at supc, and as a dip in the orbital lightcurve 
around the same phase, as shown in Fig. 19 (left panel). 
For Zo = 10^^ cm, gamma-ray absorption deforms the ob- 
served SED and the orbital lightcurve strongly, producing a 
maximum in the lightcurve around infc, = 0.72, although 
orbital changes of 9ic shift the maximum of the emission to 
slightly earlier phases. For larger Zq, e.g. Zg = 5 ■ 10^^ cm 
(Fig. 18 right panels), the effect of gamma-ray absorption 
becomes less relevant, thus only the IC angular dependence 
determines the shape of the orbital lightcurve yielding a pro- 



nounced dip around infc (see Fig. 19 right panel). In this 
case, the contribution of the counter-jet is significant around 
infc. 

We compute the gamma-ray SED and orbital 
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Figure 20. Gamma-ray SED for supc (bottom panel) and infc 
(top panel) for Vadv = 10^" cm/s. The remaining parameter 
values are i = 25°, Zq = 2 ■ lO^^ cm, r] = 10, and B = 
0.05(Z/Zo)~^ G. The jet and the counter-jet components, plus 
the summation of both, are shown. The absorbed components 
are presented with solid lines, and the production components 
are shown with dotted lines. 



lightcurves for an order of magnitude faster advection veloc- 
ity, i.e. Vadv ~ 10^" cm/s. The other parameters are fixed to 
= 2 • 10^^ cm,B = 0.05(^/Zo)"^ G, = 10, and i = 25°. 
The results are presented in Figs. [20] and |21[ To show the 
impact of Vadv, in Fig. |22]we show the ratios fast flow case 
(Vadv = 10^" cm/s) to slow flow case (Vadv = 10® cm/s) for 
the SED (left panel), and orbital lightcurves at different en- 
ergies (S^ = 0.1, 1 and 10 TeV; right panel). The SED are 
averaged over two orbital phase intervals: 0.45 < (j> < 0.9; 
(j} < 0.45 and (j) > 0.9. For the lightcurves, we focus on the 
jet radiation because the electron propagation has a weak 
impact on the counter-jet radiation. As discussed in Sec. |2.2[ 
propagation allows 1 — 10 TeV electrons to reach regions for 
which T.y.y becomes quite small. This leads to a hardening of 
the photon spectrum around supc (AF ~ 0.2). Different pat- 
terns of variability are seen for different energy bands. We 
note that for comparison, the three curves are normalised 
to the same level. At higher energies the emission becomes 



> 

CD 



■d 




0.2 



0.4 0.6 O.i 

Orbital phase 



Figure 21. Gamma-ray orbital lightcurve for the same case and 
parameter values as in Fig. |20| The jet and the counter-jet com- 
ponents, plus the summation of both, are shown. The absorbed 
components are presented with solid lines, and the production 
components are shown with dotted lines. 




Figure 22. In the left panel, the ratio fast flow case (Vadv = 
10^" cm/s) to slow flow case (Vadv = lO" cm/s) for the SED 
averaged over two orbital phase intervals: 0.45 < < 0.9 (solid 
line -INFC-); cf) < 0.45 and cf) > 0.9 (dashed line -SUPC-). We 
show added both the jet and the counter-jet component. In the 
right panel, the ratio for the same two Vadv values for the orbital 
lightcurves at three different energies (0.1 -solid line-, 1 -dotted 
line- and 10 TeV -dashed line-). The parameter values are the 
same as in Fig. |20| 



less variable because the angular dependence of the IC and 
pair-production cross-sections is weaker than at lower ener- 
gies. 



4.2 Explaining the phase averaged energy spectra 
of LS 5039 

Given the poor time resolution of the experimental data con- 
cerning the energy spectrum in LS 5039 ( i.e. high statistic 
spectra can be only given for two wide orbital phase inter- 
vals), we do not intend to make a fit to the HESS data. 
Nevertheless, below we show that, with a certain choice of 
key model parameters, we can explain the orbital phase av- 
eraged spectra reported by HESS. 

In Fig. |23[ we present the HESS spectra for two orbital 
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phase intervals (0.45 < < 0.9; < 0.45 and </> > 0.9) 
together with the gamma-ray SED calculated for the follow- 
ing parameters: T4dv = 10^° cm/s, ri = 35, So = 0.1 G, 
i — 25°, Zi) = 1.2 • lO^'^ cm, and an electron acceleration 
power Linj = 10'^^ erg/s. As is seen in Fig. 



INFC 



23 



the calculated 

SED agree reasonably with the data reported by HESS for 
both phase intervals. Our model shows changes in the flux 
by a factor ~ 2 — 3 along the orbit, similar to the observed 
ones, and the emission tends to peak around infc. Never- 
theless, we note that the lightcurves presented by HESS are 
derived from fits to the specific flux at 1 TeV. Given the un- 
certainty of the real observed spectrum, it may well be that 
the HESS fit lightcurve would be under- or overestimated 
at that energy, thus we do not carry out a comparison be- 
tween the observed and calculated lightcurves. We note that 
some other sets of parameter values could also yield a rea- 
sonable agreement with HESS data. Unfortunately, because 
of the orbital phase averaging, HESS data do not allow us to 
fix robustly the parameter space. This could be done with 
time resolved energy spectra obtained within narrow orbital 
phase intervals {/S.cf> <^ 0.1). For similar reasons, we do not 
show in Fig.[23]the data points of the EGRET source associ- 
ated to LS 5039 ( [Paredes et al.|2000 l because the data were 
taken in different epochs and do not contain orbital phase 
information. Nevertheless, it is interesting to note that the 
observed (time averaged) ratio ~GeV (EGRET) to ~TeV 
(HESS) fluxes, of about 10, is close to the model results. 



4.3 Accelerator deep inside the system 

In this section we explore the case when the accelerator is 
deep inside the system, i.e. Zq — > 0. As discussed in Sec. |2.1[ 
only an extreme accelerator with »7 < 10 can provide a 
gamma-ray spectrum extending to > 10 TeV. In Fig. |24[ we 
show the SED calculated averaging over two orbital phase 
intervals (0.45 < < 0.9; 4> < 0.45 and (/> > 0.9), taking 
i = 15, 55°, and Bq = 0.1 G. This case does not depend on 
the advection velocity. As seen in the figure, we note that 
the spectrum around supc, i.e. (p < 0.45 and (p > 0.9, can 
be barely reproduced for any inclination angle, whereas cal- 
culations for the phases around infc, 0.45 < 4> < 0.9, are 
in reasonable agreement with HESS data, in particular for 
i = 55°. 

A deep accelerator implies very large optical depths 
(r-y-y ^ 1) for photons produced all along the orbit. This 
leads to an amount of absorbed energy much larger than 
the one observed in VHE gamma-rays. However, this energy 
does not disappear, but it is released at other wavelengths. 
In particular, if the magnetic field in the system is small 



enough, ^ 10 G (see Sec. 2.4 1, absorption of gamma-rays 



will initiate an electromagnetic cascade (Aharonian et al. 
|2006b| |Bednarek|2007[ ) in the stellar radiation field. In case 
of a fully developed cascade (r ^ 1), the SED is charac- 
terised by a standard shape with a maximum at the thresh- 
old i?th ~ rri^c^/kT ~ 10 GeV, a sharp drop above Eth_, 
and a flat part at E-, > lO-Eth. Under an almost monoener- 
getic, e.g. black-body type, radiation field, the optical depth 
decreases to larger energies, thus the source becomes opti- 
cally thin and the energy spectrum at the highest energies 
is dominated by the unabsorbed intrinsic component. All 
these features can be seen in Fig. |25[ where the cascade 
SED has been calculated for the supc (non-averaged) case. 
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Figure 23. Comparison between model results and HESS ob- 
servations for both orbital phases intervals: ~infc (upper panel) 
and ~supc (lower panel). The following parameters are adopted: 
i = 25°; Zo = 1.2 ■ 10^^ cm; rj = 35; B = 0.1(Z/Zo)-^ 
Vadv = lO'^" cm/s. We show added both the jet and the counter- 
jet component. 



In LS 5039, the cascade development yields an order of mag- 
nitude or larger GeV flux compared to the TeV energy flux. 
Whereas the flat SED of the cascade radiation could explain 
the HESS spectral points in the infc phase interval, it hardly 
matches the steeper spectrum for the supc phase interval. 

If the magnetic fleld in LS 5039 exceeds 10 G, the cas- 
cade will be effectively suppressed because of dominant syn- 
chrotron cooling of the secondary electrons created every- 
where in the system. In such a case, the main fraction of 
nonthermal energy will be released in the X-ray/soft gamma- 
ray domain. This is shown in Fig. |26| in which the syn- 
chrotron radiation is shown for two values of the system av- 
erage magnetic fleld: 10 and 100 G. We note that the fluxes 
shown in Fig. [26] are calculated for an acceleration power of 
primary electrons Linj = 10'^^ erg/s, required to explain the 



TeV fluxes reported by HESS (see Fig. 241. For a system av- 
erage magnetic fleld of ~ 10 G, the synchrotron radiation of 
secondary electrons peaks around 10 keV, with a flux which 
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Figure 24. The case of a deep accelerator. Gamma-ray SED av- 
eraged over the ~supc (bottom panels) and the ~infc (top panels) 
phase intervals. Two inclination angles are adopted: i = 15° (left 
panels) and 55° (right panels). The remaining parameter values 
are Zo = 0, Vadv = 10^ cm/s, »7 = 1, and B = 0.1(Z/Zo)"^ G. 
Here, only one component (close to the compact object) is com- 
puted. 



10" 



Intrinsic 




100 



Figure 25. Gamma-ray SED obtained from a Monte Carlo calcu- 
lation of an electromagnetic cascade for the emission right at supc 
(non-averaged) in the deep accelerator case (see Fig. |24[ l. Here, 
only one component (close to the compact object) is computed. 



is close to the observed fluxes ( Bosch- Ramon et al. 20071. 
In the case of a stronger magnetic field (e.g. 100 G), the 
synchrotron emission peaks in the soft gamma-ray range, 
being also below the reported fluxes ( Strong et al.|200T l. 



5 SUMMARY 

In this paper we present the results of a detailed numerical 
study of the spectral and temporal properties of the high en- 
ergy gamma-ray emission from the binary system LS 5039 
adopting an IC scenario. The calculations have been con- 
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Figure 26. Synchrotron radiation of secondary electrons for two 
magnetic fields: iJ = 10 (dashed line), 100 G (solid line). 



ducted in the framework of a specific model which assumes 
that gamma-rays are produced in a jet-like structure perpen- 
dicular to the orbital plane. However, the obtained results 
and conclusions concerning the electron accelerator and the 
gamma-ray emitter have broader implications. The sheer 
fact that we see a modulated gamma-ray signal with a 
3.9 day period, with a hard energy spectrum extending to 
10 TeV and beyond, implies that: (a) the gamma-ray emitter 
is located close to the binary system, i.e. at distances com- 
parable to RoTh\ and (b) we deal with an efficient accelerator 
with acceleration rate not far from the theoretically max- 
imum possible rate, i.e. rj < 100 {r/ < 10 for an accelerator 
located deep inside the binary system). Another interesting 
model-independent conclusion is that if the gamma rays 
have an inverse Compton origin, then we can restrict the 
magnetic field in the gamma-ray emitter within a quite nar- 
row range. The fact that the IC scattering which produces 
the TeV photons occurs deep in the KN regime implies that 
the observed spectrum with photon index T ~ 2 requires 
a very hard electron energy distribution with power-law in- 
dex ~ 1. Such a hard spectrum cannot be formed even in 
the case of a monoenergetic (e.g. Maxwellian type) injection 
electron spectrum unless they cool through the IC scatter- 
ing in the KN regim«[^ (see Eq.l]?])). On the other hand the 
magnetic field cannot be much smaller than 0.1 G in order 
to provide an adequate acceleration rate. Thus, the accel- 
erator emitter magnetic field should be confined in a band 
around 0.1 G. 

Because of the close location of the gamma-ray emit- 
ter to the companion star, the modification of the spec- 
trum of gamma-rays due to the gamma-gamma absorp- 
tion process is unavoidable. Moreover, gamma-ray absorp- 
tion strongly depends on the orbital phase. Because of the 
slightly elliptical orbit, the phase-dependence of the pho- 
ton density is moderate. However, the orbital change of 
the interaction angle results in a change in the threshold 
of the interaction, and thus leads to strong modulation of 

^ Synchrotron/Thomson IC cooling would yield a oc Eg^ electron 
energy distribution. 
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the optical depth. Interestingly, the energy-dependent mod- 
ification of the gamma-ray energy spectrum produced by 
gamma-gamma absorption cannot explain the data reported 
by HESS for two phase intervals. In particular, whereas one 
should expect the strongest variation of the gamma-ray flux 
at ~ 100 GeV; HESS results show just the opposite. Nev- 
ertheless, the anisotropy of the radiation field has a strong 
impact not only on the gamma-gamma absorption, but also 
on the spectral characteristics of the IC scattering. We show 
here that the combination of these two effects, both related 
to the anisotropic character of the interaction of the stellar 
radiation field with relativistic electrons (IC) and gamma- 
rays (pair production), can satisfactorily explain the HESS 
observations reported for the supc and infc phase intervals. 
In addition, we note that the phase-dependence of the elec- 
tron maximum energy and the electron advection along the 
jet could play an important role in the formation of the 
gamma-ray spectra. 

The absorption of gamma-rays does not imply that this 
energy disappears. The secondary electrons interacting with 
the stellar radiation and ambient magnetic fields produce 
secondary radiation. There are two well defined possibil- 
ities; either the energy density of the ambient medium is 
dominated by the stellar radiation field or by the magnetic 
field. In the former case, pair production and further IC 
scattering initiate an electromagnetic cascade. In the latter 
case, synchrotron radiation is produced. As noted above, 
a typical ambient magnetic field determining the boundary 
between these two regimes is several Gauss. When an effi- 
cient cascade occurs, the source becomes more transparent 
at very high energies, although the energy is mainly re- 
leased around the threshold energy ~ 10 GeV. In the case 
of the suppression of the cascade through synchrotron cool- 
ing, the source is less transparent, and synchrotron radiation 
is released mostly at keV-MeV energies, with luminosities 
exceeding the VHE gamma-ray luminosity. Interestingly, a 
~ 10 G magnetic field is plausible for the surroundings of 
a massive star, so a priori one cannot exclude any of both 
possible scenarios, which in fact may be realised simultane- 
ously in different parts of the system. In such a case, both 
cascading and synchrotron radiation of secondaries would 
contribute comparably to the emerging high energy radia- 
tion. In this regard, future deep X-ray and gamma-ray obser- 
vations, together with TeV observations with more sensitive 
instruments capable of spectral measurements on short time 
scales, will provide insights into the origin of the nonthermal 
emission of this enigmatic source. 
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